guides(color=guide_legend(byrow=TRUE))
#PLOT DATA - ACID VS NO2
ggplot(df3_427, aes(x=as.numeric(NO2), y=AbsValue, group=as.factor(Acid))) +
geom_point(aes(color=Acid), size = 2, alpha=0.3) +
stat_summary(fun.data = "mean_sdl", aes(color=Acid), size = 1, width=0.002, geom="errorbar") +
stat_summary(fun.data = "mean_sdl", aes(color=Acid), size = 1, geom="line") +
stat_summary(fun.data = "mean_sdl", aes(color=Acid), size = 4, geom="point") +
labs(x = expression("NO"[2] ~ "(%)"),
y = "Absorbance 427 nm (a.u.)") +
scale_y_continuous(breaks = seq(0, 0.4, by = 0.05)) +
scale_x_continuous(labels = scales::percent, limits = c(0,0.08)) +
theme_classic() +
theme(aspect.ratio=1,
axis.title=element_text(size=18),
axis.text=element_text(colour = 1, size = 14),
legend.title=element_blank(),
legend.justification=c(0,1),
legend.position=c(0.8,0.9),
legend.text=element_text(size = 18),
panel.border = element_rect(colour = "black", fill=NA)) +
scale_color_discrete(labels=c(expression("HBF"[4]), "HCl")) +
guides(color=guide_legend(byrow=TRUE))
#PLOT DATA - ACID VS NO2
ggplot(df3_427, aes(x=as.numeric(NO2), y=AbsValue, group=as.factor(Acid))) +
geom_point(aes(color=Acid), size = 2, alpha=0.3) +
stat_summary(fun.data = "mean_sdl", aes(color=Acid), size = 1, width=0.002, geom="errorbar") +
stat_summary(fun.data = "mean_sdl", aes(color=Acid), size = 1, geom="line") +
stat_summary(fun.data = "mean_sdl", aes(color=Acid), size = 4, geom="point") +
labs(x = expression("NO"[2] ~ "(%)"),
y = "Absorbance 427 nm (a.u.)") +
scale_y_continuous(breaks = seq(0, 0.4, by = 0.05)) +
scale_x_continuous(labels = scales::percent, limits = c(0,0.08)) +
theme_classic() +
theme(aspect.ratio=1,
axis.title=element_text(size=18),
axis.text=element_text(colour = 1, size = 14),
legend.title=element_blank(),
legend.justification=c(0,1),
legend.position=c(0.8,0.9),
legend.text=element_text(size = 18),
panel.border = element_rect(colour = "black", fill=NA)) +
scale_color_discrete(labels=c(expression("HBF"[4]), "HCl")) +
guides(color=guide_legend(byrow=TRUE, nrow = 1))
#PLOT DATA - ACID VS NO2
ggplot(df3_427, aes(x=as.numeric(NO2), y=AbsValue, group=as.factor(Acid))) +
geom_point(aes(color=Acid), size = 2, alpha=0.3) +
stat_summary(fun.data = "mean_sdl", aes(color=Acid), size = 1, width=0.002, geom="errorbar") +
stat_summary(fun.data = "mean_sdl", aes(color=Acid), size = 1, geom="line") +
stat_summary(fun.data = "mean_sdl", aes(color=Acid), size = 4, geom="point") +
labs(x = expression("NO"[2] ~ "(%)"),
y = "Absorbance 427 nm (a.u.)") +
scale_y_continuous(breaks = seq(0, 0.4, by = 0.05)) +
scale_x_continuous(labels = scales::percent, limits = c(0,0.06)) +
theme_classic() +
theme(aspect.ratio=1,
axis.title=element_text(size=18),
axis.text=element_text(colour = 1, size = 14),
legend.title=element_blank(),
legend.justification=c(0,1),
legend.position=c(0.05,0.9),
legend.text=element_text(size = 18),
panel.border = element_rect(colour = "black", fill=NA)) +
scale_color_discrete(labels=c(expression("HBF"[4]), "HCl")) +
guides(color=guide_legend(byrow=TRUE, nrow = 1))
#PLOT DATA - ACID VS NO2
ggplot(df3_427, aes(x=as.numeric(NO2), y=AbsValue, group=as.factor(Acid))) +
geom_point(aes(color=Acid), size = 2, alpha=0.3) +
stat_summary(fun.data = "mean_sdl", aes(color=Acid), size = 1, width=0.002, geom="errorbar") +
stat_summary(fun.data = "mean_sdl", aes(color=Acid), size = 1, geom="line") +
stat_summary(fun.data = "mean_sdl", aes(color=Acid), size = 4, geom="point") +
labs(x = expression("NO"[2] ~ "(%)"),
y = "Absorbance 427 nm (a.u.)") +
scale_y_continuous(breaks = seq(0, 0.4, by = 0.05)) +
scale_x_continuous(labels = scales::percent, limits = c(0,0.06)) +
theme_classic() +
theme(aspect.ratio=1,
axis.title=element_text(size=18),
axis.text=element_text(colour = 1, size = 14),
legend.title=element_blank(),
legend.justification=c(0,1),
legend.position=c(0.05,0.95),
legend.text=element_text(size = 18),
panel.border = element_rect(colour = "black", fill=NA)) +
scale_color_discrete(labels=c(expression("HBF"[4]), "HCl")) +
guides(color=guide_legend(byrow=TRUE, nrow = 1))
#PLOT DATA - ACID VS NO2
ggplot(df3_427, aes(x=as.numeric(NO2), y=AbsValue, group=as.factor(Acid))) +
geom_point(aes(color=Acid), size = 2, alpha=0.3) +
stat_summary(fun.data = "mean_sdl", aes(color=Acid), size = 1, width=0.002, geom="errorbar") +
stat_summary(fun.data = "mean_sdl", aes(color=Acid), size = 1, geom="line") +
stat_summary(fun.data = "mean_sdl", aes(color=Acid), size = 4, geom="point") +
labs(x = expression("NO"[2] ~ "(%)"),
y = "Absorbance 427 nm (a.u.)") +
scale_y_continuous(breaks = seq(0, 0.4, by = 0.05)) +
scale_x_continuous(labels = scales::percent, limits = c(0,0.06)) +
theme_classic() +
theme(aspect.ratio=1,
axis.title=element_text(size=18),
axis.text=element_text(colour = 1, size = 14),
legend.title=element_blank(),
legend.justification=c(0,1),
legend.position=c(0.05,0.97),
legend.text=element_text(size = 18),
panel.border = element_rect(colour = "black", fill=NA)) +
scale_color_discrete(labels=c(expression("HBF"[4]), "HCl")) +
guides(color=guide_legend(byrow=TRUE, nrow = 1))
#PLOT DATA - ACID VS NO2
ggplot(df3_427, aes(x=as.numeric(NO2), y=AbsValue, group=as.factor(Acid))) +
geom_point(aes(color=Acid), size = 2, alpha=0.3) +
stat_summary(fun.data = "mean_sdl", aes(color=Acid), size = 1, width=0.002, geom="errorbar") +
stat_summary(fun.data = "mean_sdl", aes(color=Acid), size = 1, geom="line") +
stat_summary(fun.data = "mean_sdl", aes(color=Acid), size = 4, geom="point") +
labs(x = expression("NO"[2] ~ "(%)"),
y = "Absorbance 427 nm (a.u.)") +
scale_y_continuous(breaks = seq(0, 0.4, by = 0.05)) +
scale_x_continuous(labels = scales::percent, limits = c(0,0.06)) +
theme_classic() +
theme(aspect.ratio=1,
axis.title=element_text(size=18),
axis.text=element_text(colour = 1, size = 14),
legend.title=element_blank(),
legend.justification=c(0,1),
legend.position=c(0.05,0.97),
legend.text=element_text(size = 18),
panel.border = element_rect(colour = "black", fill=NA)) +
scale_color_discrete(labels=c(expression("HBF"[4]), "HCl")) +
guides(color=guide_legend(byrow=TRUE, nrow = 1))
#PRINT/SAVE IMAGES, CHANGE FILE NAMES AS APPROPRIATE
ggsave("FactorScreen_NO2vsAcid.pdf", height=13, width=13, units="cm")
#PLOT DATA - ACID VS NO2
ggplot(df3_427, aes(x=as.numeric(NO2), y=AbsValue, group=as.factor(Acid))) +
geom_point(aes(color=Acid), size = 2, alpha=0.3) +
stat_summary(fun.data = "mean_sdl", aes(color=Acid), size = 1, width=0.002, geom="errorbar") +
stat_summary(fun.data = "mean_sdl", aes(color=Acid), size = 1, geom="line") +
stat_summary(fun.data = "mean_sdl", aes(color=Acid), size = 4, geom="point") +
labs(x = expression("NO"[2] ~ "(%)"),
y = "Absorbance 427 nm (a.u.)") +
scale_y_continuous(breaks = seq(0, 0.4, by = 0.05)) +
scale_x_continuous(labels = scales::percent, limits = c(0,0.06)) +
theme_classic() +
theme(aspect.ratio=1,
axis.title=element_text(size=18),
axis.text=element_text(colour = 1, size = 14),
legend.title=element_blank(),
legend.justification=c(0,1),
legend.position=c(0.05,0.97),
legend.text=element_text(size = 18),
panel.border = element_rect(colour = "black", fill=NA))
#PLOT DATA - ACID VS AMTA
ggplot(df3_427, aes(x=AMTA, y=AbsValue, group=as.factor(Acid))) +
geom_point(aes(color=Acid), size = 2, alpha=0.3) +
stat_summary(fun.data = "mean_sdl", aes(color=Acid), size = 1, width=0.2, geom="errorbar") +
stat_summary(fun.data = "mean_sdl", aes(color=Acid), size = 1, geom="line") +
stat_summary(fun.data = "mean_sdl", aes(color=Acid), size = 4, geom="point") +
labs(x = "[AMTA] (mM)",
y = "Absorbance 427 nm (a.u.)") +
scale_y_continuous(breaks = seq(0, 0.4, by = 0.05)) +
theme_classic() +
theme(aspect.ratio=1,
axis.title=element_text(size=18),
axis.text=element_text(colour = 1, size = 14),
legend.title=element_blank(),
legend.justification=c(0,1),
legend.position=c(0.05,0.9),
legend.text=element_text(size = 18),
panel.border = element_rect(colour = "black", fill=NA)) +
scale_color_discrete(labels=c(expression("HBF"[4]), "HCl")) +
guides(color=guide_legend(byrow=TRUE, nrow = 1))
#PLOT DATA - ACID VS AMTA
ggplot(df3_427, aes(x=AMTA, y=AbsValue, group=as.factor(Acid))) +
geom_point(aes(color=Acid), size = 2, alpha=0.3) +
stat_summary(fun.data = "mean_sdl", aes(color=Acid), size = 1, width=0.2, geom="errorbar") +
stat_summary(fun.data = "mean_sdl", aes(color=Acid), size = 1, geom="line") +
stat_summary(fun.data = "mean_sdl", aes(color=Acid), size = 4, geom="point") +
labs(x = "[AMTA] (mM)",
y = "Absorbance 427 nm (a.u.)") +
scale_y_continuous(breaks = seq(0, 0.4, by = 0.05)) +
theme_classic() +
theme(aspect.ratio=1,
axis.title=element_text(size=18),
axis.text=element_text(colour = 1, size = 14),
legend.title=element_blank(),
legend.justification=c(0,1),
legend.position=c(0.05,0.97),
legend.text=element_text(size = 18),
panel.border = element_rect(colour = "black", fill=NA)) +
scale_color_discrete(labels=c(expression("HBF"[4]), "HCl")) +
guides(color=guide_legend(byrow=TRUE, nrow = 1))
#PRINT/SAVE IMAGES, CHANGE FILE NAMES AS APPROPRIATE
ggsave("FactorScreen_AMTAvsAcid.pdf", height=13, width=13, units="cm")
#PLOT DATA - ACID VS AMTA
ggplot(df3_427, aes(x=AMTA, y=AbsValue, group=as.factor(Acid))) +
geom_point(aes(color=Acid), size = 2, alpha=0.3) +
stat_summary(fun.data = "mean_sdl", aes(color=Acid), size = 1, width=0.2, geom="errorbar") +
stat_summary(fun.data = "mean_sdl", aes(color=Acid), size = 1, geom="line") +
stat_summary(fun.data = "mean_sdl", aes(color=Acid), size = 4, geom="point") +
labs(x = "[AMTA] (mM)",
y = "Absorbance 427 nm (a.u.)") +
scale_y_continuous(breaks = seq(-0.02, 0.4, by = 0.05)) +
theme_classic() +
theme(aspect.ratio=1,
axis.title=element_text(size=18),
axis.text=element_text(colour = 1, size = 14),
legend.title=element_blank(),
legend.justification=c(0,1),
legend.position=c(0.05,0.97),
legend.text=element_text(size = 18),
panel.border = element_rect(colour = "black", fill=NA)) +
scale_color_discrete(labels=c(expression("HBF"[4]), "HCl")) +
guides(color=guide_legend(byrow=TRUE, nrow = 1))
#PLOT DATA - ACID VS AMTA
ggplot(df3_427, aes(x=AMTA, y=AbsValue, group=as.factor(Acid))) +
geom_point(aes(color=Acid), size = 2, alpha=0.3) +
stat_summary(fun.data = "mean_sdl", aes(color=Acid), size = 1, width=0.2, geom="errorbar") +
stat_summary(fun.data = "mean_sdl", aes(color=Acid), size = 1, geom="line") +
stat_summary(fun.data = "mean_sdl", aes(color=Acid), size = 4, geom="point") +
labs(x = "[AMTA] (mM)",
y = "Absorbance 427 nm (a.u.)") +
scale_y_continuous(breaks = seq(0, 0.4, by = 0.05)) +
theme_classic() +
theme(aspect.ratio=1,
axis.title=element_text(size=18),
axis.text=element_text(colour = 1, size = 14),
legend.title=element_blank(),
legend.justification=c(0,1),
legend.position=c(0.05,0.97),
legend.text=element_text(size = 18),
panel.border = element_rect(colour = "black", fill=NA)) +
scale_color_discrete(labels=c(expression("HBF"[4]), "HCl")) +
guides(color=guide_legend(byrow=TRUE, nrow = 1))
#PRINT/SAVE IMAGES, CHANGE FILE NAMES AS APPROPRIATE
ggsave("FactorScreen_AMTAvsAcid.pdf", height=13, width=13, units="cm")
#PLOT DATA - ACID VS NO2
ggplot(df3_427, aes(x=as.numeric(NO2), y=AbsValue, group=as.factor(Acid))) +
geom_point(aes(color=Acid), size = 2, alpha=0.3) +
stat_summary(fun.data = "mean_sdl", aes(color=Acid), size = 1, width=0.002, geom="errorbar") +
stat_summary(fun.data = "mean_sdl", aes(color=Acid), size = 1, geom="line") +
stat_summary(fun.data = "mean_sdl", aes(color=Acid), size = 4, geom="point") +
labs(x = expression("NO"[2] ~ "(%)"),
y = "Absorbance 427 nm (a.u.)") +
scale_y_continuous(breaks = seq(0, 0.4, by = 0.05)) +
scale_x_continuous(labels = scales::percent, limits = c(0,0.06)) +
theme_classic() +
theme(aspect.ratio=1,
axis.title=element_text(size=18),
axis.text=element_text(colour = 1, size = 14),
legend.title=element_blank(),
legend.justification=c(0,1),
legend.position=c(0.05,0.97),
legend.text=element_text(size = 18),
panel.border = element_rect(colour = "black", fill=NA)) +
scale_color_discrete(labels=c(expression("HBF"[4]), "HCl")) +
guides(color=guide_legend(byrow=TRUE, nrow = 1))
#PRINT/SAVE IMAGES, CHANGE FILE NAMES AS APPROPRIATE
ggsave("FactorScreen_NO2vsAcid.pdf", height=13, width=13, units="cm")
#PLOT DATA - AMTA VS NO2
ggplot(df3_427, aes(x=as.numeric(NO2), y=AbsValue, group=as.factor(AMTA))) +
geom_point(aes(color=Acid), size = 2, alpha=0.3) +
stat_summary(fun.data = "mean_sdl", aes(color=AMTA), size = 1, width=0.002, geom="errorbar") +
stat_summary(fun.data = "mean_sdl", aes(color=AMTA), size = 1, geom="line") +
stat_summary(fun.data = "mean_sdl", aes(color=AMTA), size = 4, geom="point") +
labs(x = expression("NO"[2] ~ "(%)"),
y = "Absorbance 427 nm (a.u.)") +
scale_y_continuous(breaks = seq(0, 0.4, by = 0.05)) +
scale_x_continuous(labels = scales::percent, limits = c(0,0.06)) +
theme_classic() +
theme(aspect.ratio=1,
axis.title=element_text(size=18),
axis.text=element_text(colour = 1, size = 14),
legend.title=element_blank(),
legend.justification=c(0,1),
legend.position=c(0.05,0.97),
legend.text=element_text(size = 18),
panel.border = element_rect(colour = "black", fill=NA)) +
scale_color_discrete(labels=c(expression("HBF"[4]), "HCl")) +
guides(color=guide_legend(byrow=TRUE, nrow = 1))
#PLOT DATA - AMTA VS NO2
ggplot(df3_427, aes(x=as.numeric(NO2), y=AbsValue, group=as.factor(AMTA))) +
geom_point(aes(color=Acid), size = 2, alpha=0.3) +
stat_summary(fun.data = "mean_sdl", aes(color=AMTA), size = 1, width=0.002, geom="errorbar") +
stat_summary(fun.data = "mean_sdl", aes(color=AMTA), size = 1, geom="line") +
stat_summary(fun.data = "mean_sdl", aes(color=AMTA), size = 4, geom="point") +
labs(x = expression("NO"[2] ~ "(%)"),
y = "Absorbance 427 nm (a.u.)") +
scale_y_continuous(breaks = seq(0, 0.4, by = 0.05)) +
scale_x_continuous(labels = scales::percent, limits = c(0,0.06)) +
theme_classic() +
theme(aspect.ratio=1,
axis.title=element_text(size=18),
axis.text=element_text(colour = 1, size = 14),
legend.title=element_blank(),
legend.justification=c(0,1),
legend.position=c(0.05,0.97),
legend.text=element_text(size = 18),
panel.border = element_rect(colour = "black", fill=NA)) +
guides(color=guide_legend(byrow=TRUE, nrow = 1))
#PLOT DATA - AMTA VS NO2
ggplot(df3_427, aes(x=as.numeric(NO2), y=AbsValue, group=as.factor(AMTA))) +
geom_point(aes(color=AMTA), size = 2, alpha=0.3) +
stat_summary(fun.data = "mean_sdl", aes(color=AMTA), size = 1, width=0.002, geom="errorbar") +
stat_summary(fun.data = "mean_sdl", aes(color=AMTA), size = 1, geom="line") +
stat_summary(fun.data = "mean_sdl", aes(color=AMTA), size = 4, geom="point") +
labs(x = expression("NO"[2] ~ "(%)"),
y = "Absorbance 427 nm (a.u.)") +
scale_y_continuous(breaks = seq(0, 0.4, by = 0.05)) +
scale_x_continuous(labels = scales::percent, limits = c(0,0.06)) +
theme_classic() +
theme(aspect.ratio=1,
axis.title=element_text(size=18),
axis.text=element_text(colour = 1, size = 14),
legend.title=element_blank(),
legend.justification=c(0,1),
legend.position=c(0.05,0.97),
legend.text=element_text(size = 18),
panel.border = element_rect(colour = "black", fill=NA)) +
guides(color=guide_legend(byrow=TRUE, nrow = 1))
#PLOT DATA - AMTA VS NO2
ggplot(df3_427, aes(x=as.numeric(NO2), y=AbsValue, group=as.factor(AMTA))) +
geom_point(aes(color=AMTA), size = 2, alpha=0.3) +
stat_summary(fun.data = "mean_sdl", aes(color=AMTA), size = 1, width=0.002, geom="errorbar") +
stat_summary(fun.data = "mean_sdl", aes(color=AMTA), size = 1, geom="line") +
stat_summary(fun.data = "mean_sdl", aes(color=AMTA), size = 4, geom="point") +
labs(x = expression("NO"[2] ~ "(%)"),
y = "Absorbance 427 nm (a.u.)") +
scale_y_continuous(breaks = seq(0, 0.4, by = 0.05)) +
scale_x_continuous(labels = scales::percent, limits = c(0,0.06)) +
theme_classic() +
theme(aspect.ratio=1,
axis.title=element_text(size=18),
axis.text=element_text(colour = 1, size = 14),
legend.title=element_blank(),
legend.justification=c(0,1),
legend.position=c(0.05,0.97),
legend.text=element_text(size = 18),
panel.border = element_rect(colour = "black", fill=NA)) +
scale_color_discrete(labels=c("0.1 mM AMTA", "1 mM AMTA")) +
guides(color=guide_legend(byrow=TRUE, nrow = 1))
#PRINT/SAVE IMAGES, CHANGE FILE NAMES AS APPROPRIATE
ggsave("FactorScreen_NO2vsAMTA.pdf", height=13, width=13, units="cm")
#PLOT DATA - AMTA VS NO2
ggplot(df3_427, aes(x=as.numeric(NO2), y=AbsValue, group=as.factor(AMTA))) +
geom_point(aes(color=AMTA), size = 2, alpha=0.3) +
stat_summary(fun.data = "mean_sdl", aes(color=AMTA), size = 1, width=0.002, geom="errorbar") +
stat_summary(fun.data = "mean_sdl", aes(color=AMTA), size = 1, geom="line") +
stat_summary(fun.data = "mean_sdl", aes(color=AMTA), size = 4, geom="point") +
labs(x = expression("NO"[2] ~ "(%)"),
y = "Absorbance 427 nm (a.u.)") +
scale_y_continuous(breaks = seq(0, 0.4, by = 0.05), limits = c(-0.05,0.3)) +
scale_x_continuous(labels = scales::percent, limits = c(0,0.06)) +
theme_classic() +
theme(aspect.ratio=1,
axis.title=element_text(size=18),
axis.text=element_text(colour = 1, size = 14),
legend.title=element_blank(),
legend.justification=c(0,1),
legend.position=c(0.05,0.97),
legend.text=element_text(size = 18),
panel.border = element_rect(colour = "black", fill=NA)) +
scale_color_discrete(labels=c("0.1 mM AMTA", "1 mM AMTA")) +
guides(color=guide_legend(byrow=TRUE, nrow = 1))
#PLOT DATA - AMTA VS NO2
ggplot(df3_427, aes(x=as.numeric(NO2), y=AbsValue, group=as.factor(AMTA))) +
geom_point(aes(color=AMTA), size = 2, alpha=0.3) +
stat_summary(fun.data = "mean_sdl", aes(color=AMTA), size = 1, width=0.002, geom="errorbar") +
stat_summary(fun.data = "mean_sdl", aes(color=AMTA), size = 1, geom="line") +
stat_summary(fun.data = "mean_sdl", aes(color=AMTA), size = 4, geom="point") +
labs(x = expression("NO"[2] ~ "(%)"),
y = "Absorbance 427 nm (a.u.)") +
scale_y_continuous(breaks = seq(0, 0.4, by = 0.05), limits = c(-0.02,0.3)) +
scale_x_continuous(labels = scales::percent, limits = c(0,0.06)) +
theme_classic() +
theme(aspect.ratio=1,
axis.title=element_text(size=18),
axis.text=element_text(colour = 1, size = 14),
legend.title=element_blank(),
legend.justification=c(0,1),
legend.position=c(0.05,0.97),
legend.text=element_text(size = 18),
panel.border = element_rect(colour = "black", fill=NA)) +
scale_color_discrete(labels=c("0.1 mM AMTA", "1 mM AMTA")) +
guides(color=guide_legend(byrow=TRUE, nrow = 1))
#PLOT DATA - AMTA VS NO2
ggplot(df3_427, aes(x=as.numeric(NO2), y=AbsValue, group=as.factor(AMTA))) +
geom_point(aes(color=AMTA), size = 2, alpha=0.3) +
stat_summary(fun.data = "mean_sdl", aes(color=AMTA), size = 1, width=0.002, geom="errorbar") +
stat_summary(fun.data = "mean_sdl", aes(color=AMTA), size = 1, geom="line") +
stat_summary(fun.data = "mean_sdl", aes(color=AMTA), size = 4, geom="point") +
labs(x = expression("NO"[2] ~ "(%)"),
y = "Absorbance 427 nm (a.u.)") +
scale_y_continuous(breaks = seq(0, 0.4, by = 0.05), limits = c(-0.02,0.28)) +
scale_x_continuous(labels = scales::percent, limits = c(0,0.06)) +
theme_classic() +
theme(aspect.ratio=1,
axis.title=element_text(size=18),
axis.text=element_text(colour = 1, size = 14),
legend.title=element_blank(),
legend.justification=c(0,1),
legend.position=c(0.05,0.97),
legend.text=element_text(size = 18),
panel.border = element_rect(colour = "black", fill=NA)) +
scale_color_discrete(labels=c("0.1 mM AMTA", "1 mM AMTA")) +
guides(color=guide_legend(byrow=TRUE, nrow = 1))
#PLOT DATA - ACID VS NO2
ggplot(df3_427, aes(x=as.numeric(NO2), y=AbsValue, group=as.factor(Acid))) +
geom_point(aes(color=Acid), size = 2, alpha=0.3) +
stat_summary(fun.data = "mean_sdl", aes(color=Acid), size = 1, width=0.002, geom="errorbar") +
stat_summary(fun.data = "mean_sdl", aes(color=Acid), size = 1, geom="line") +
stat_summary(fun.data = "mean_sdl", aes(color=Acid), size = 4, geom="point") +
labs(x = expression("NO"[2] ~ "(%)"),
y = "Absorbance 427 nm (a.u.)") +
scale_y_continuous(breaks = seq(0, 0.4, by = 0.05), limits = c(-0.02,0.28)) +
scale_x_continuous(labels = scales::percent, limits = c(0,0.06)) +
theme_classic() +
theme(aspect.ratio=1,
axis.title=element_text(size=18),
axis.text=element_text(colour = 1, size = 14),
legend.title=element_blank(),
legend.justification=c(0,1),
legend.position=c(0.05,0.97),
legend.text=element_text(size = 18),
panel.border = element_rect(colour = "black", fill=NA)) +
scale_color_discrete(labels=c(expression("HBF"[4]), "HCl")) +
guides(color=guide_legend(byrow=TRUE, nrow = 1))
#PLOT DATA - AMTA VS NO2
ggplot(df3_427, aes(x=as.numeric(NO2), y=AbsValue, group=as.factor(AMTA))) +
geom_point(aes(color=AMTA), size = 2, alpha=0.3) +
stat_summary(fun.data = "mean_sdl", aes(color=AMTA), size = 1, width=0.002, geom="errorbar") +
stat_summary(fun.data = "mean_sdl", aes(color=AMTA), size = 1, geom="line") +
stat_summary(fun.data = "mean_sdl", aes(color=AMTA), size = 4, geom="point") +
labs(x = expression("NO"[2] ~ "(%)"),
y = "Absorbance 427 nm (a.u.)") +
scale_y_continuous(breaks = seq(0, 0.4, by = 0.05), limits = c(-0.02,0.28)) +
scale_x_continuous(labels = scales::percent, limits = c(0,0.06)) +
theme_classic() +
theme(aspect.ratio=1,
axis.title=element_text(size=18),
axis.text=element_text(colour = 1, size = 14),
legend.title=element_blank(),
legend.justification=c(0,1),
legend.position=c(0.05,0.97),
legend.text=element_text(size = 18),
panel.border = element_rect(colour = "black", fill=NA)) +
scale_color_discrete(labels=c("0.1 mM AMTA", "1 mM AMTA")) +
guides(color=guide_legend(byrow=TRUE, nrow = 1))
#PLOT DATA - ACID VS AMTA
ggplot(df3_427, aes(x=AMTA, y=AbsValue, group=as.factor(Acid))) +
geom_point(aes(color=Acid), size = 2, alpha=0.3) +
stat_summary(fun.data = "mean_sdl", aes(color=Acid), size = 1, width=0.2, geom="errorbar") +
stat_summary(fun.data = "mean_sdl", aes(color=Acid), size = 1, geom="line") +
stat_summary(fun.data = "mean_sdl", aes(color=Acid), size = 4, geom="point") +
labs(x = "[AMTA] (mM)",
y = "Absorbance 427 nm (a.u.)") +
scale_y_continuous(breaks = seq(0, 0.4, by = 0.05), limits = c(-0.02,0.28)) +
theme_classic() +
theme(aspect.ratio=1,
axis.title=element_text(size=18),
axis.text=element_text(colour = 1, size = 14),
legend.title=element_blank(),
legend.justification=c(0,1),
legend.position=c(0.05,0.97),
legend.text=element_text(size = 18),
panel.border = element_rect(colour = "black", fill=NA)) +
scale_color_discrete(labels=c(expression("HBF"[4]), "HCl")) +
guides(color=guide_legend(byrow=TRUE, nrow = 1))
#PRINT/SAVE IMAGES, CHANGE FILE NAMES AS APPROPRIATE
ggsave("FactorScreen_AMTAvsAcid.pdf", height=13, width=13, units="cm")
#PLOT DATA - ACID VS NO2
ggplot(df3_427, aes(x=as.numeric(NO2), y=AbsValue, group=as.factor(Acid))) +
geom_point(aes(color=Acid), size = 2, alpha=0.3) +
stat_summary(fun.data = "mean_sdl", aes(color=Acid), size = 1, width=0.002, geom="errorbar") +
stat_summary(fun.data = "mean_sdl", aes(color=Acid), size = 1, geom="line") +
stat_summary(fun.data = "mean_sdl", aes(color=Acid), size = 4, geom="point") +
labs(x = expression("NO"[2] ~ "(%)"),
y = "Absorbance 427 nm (a.u.)") +
scale_y_continuous(breaks = seq(0, 0.4, by = 0.05), limits = c(-0.02,0.28)) +
scale_x_continuous(labels = scales::percent, limits = c(0,0.06)) +
theme_classic() +
theme(aspect.ratio=1,
axis.title=element_text(size=18),
axis.text=element_text(colour = 1, size = 14),
legend.title=element_blank(),
legend.justification=c(0,1),
legend.position=c(0.05,0.97),
legend.text=element_text(size = 18),
panel.border = element_rect(colour = "black", fill=NA)) +
scale_color_discrete(labels=c(expression("HBF"[4]), "HCl")) +
guides(color=guide_legend(byrow=TRUE, nrow = 1))
#PRINT/SAVE IMAGES, CHANGE FILE NAMES AS APPROPRIATE
ggsave("FactorScreen_NO2vsAcid.pdf", height=13, width=13, units="cm")
#PLOT DATA - AMTA VS NO2
ggplot(df3_427, aes(x=as.numeric(NO2), y=AbsValue, group=as.factor(AMTA))) +
geom_point(aes(color=AMTA), size = 2, alpha=0.3) +
stat_summary(fun.data = "mean_sdl", aes(color=AMTA), size = 1, width=0.002, geom="errorbar") +
stat_summary(fun.data = "mean_sdl", aes(color=AMTA), size = 1, geom="line") +
stat_summary(fun.data = "mean_sdl", aes(color=AMTA), size = 4, geom="point") +
labs(x = expression("NO"[2] ~ "(%)"),
y = "Absorbance 427 nm (a.u.)") +
scale_y_continuous(breaks = seq(0, 0.4, by = 0.05), limits = c(-0.02,0.28)) +
scale_x_continuous(labels = scales::percent, limits = c(0,0.06)) +
theme_classic() +
theme(aspect.ratio=1,
axis.title=element_text(size=18),
axis.text=element_text(colour = 1, size = 14),
legend.title=element_blank(),
legend.justification=c(0,1),
legend.position=c(0.05,0.97),
legend.text=element_text(size = 18),
panel.border = element_rect(colour = "black", fill=NA)) +
scale_color_discrete(labels=c("0.1 mM AMTA", "1 mM AMTA")) +
guides(color=guide_legend(byrow=TRUE, nrow = 1))
#PRINT/SAVE IMAGES, CHANGE FILE NAMES AS APPROPRIATE
ggsave("FactorScreen_NO2vsAMTA.pdf", height=13, width=13, units="cm")
